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Topological objects of intricate structures have been found in a wide range of systems, 
like cosmology to liquid crystal^, DNA chains^, superfluid '^He (Ref. 3), quantum Hall 
magnets^'^, Bose-Einstein condensates^ etc. The topological spin texture in helical mag- 
nets, observed in recently^'^, makes a new entry into this fascinating phenomenon. The 
unusual magnetic behaviour of helimagnet MnSi, noticed in recent years^^^^, prompted the 
suspicion that the magnetic states arising in such crystals are of topological nature^'^^^. 
Experiments''^ based on the topological Hall effect confirmed such topologically nontrivial 
states as the skyrmions^^, located on a plane^^'^^ perpendicular to the applied magnetic 
field. However, the available models close to MnSi, investigating the formation of skyrmion 
states^'^'^'^^'^^, are based mostly on approximate or numerical methods. We present here a 
theoretical model for chiral magnets with competing exchange and Dzyaloshinskii-Moriya 
type interaction, which leads to an exact skyrmionic solution with integer topological charge 
A^. Such topologically stable spin states with analytic solution on a two-dimensional plane^^ 
show helical structures of partial order, without inversion and circular symmetry. These 
exact N-skyrmions, though represent higher excited states, correspond to the lowest en- 
ergy stable configuration in each topological sector and are likely to appear in MnSi under 
suitable experimental conditions. The present exact topological solitons, with explicit non- 
circular symmetry could be applicable also to other fields, where skyrmions are observed, 
especially in natural systems with less symmetries. 

Topological properties can be revealed through a mapping from a continuum space to a differentiable 
manifold. Therefore for describing topological objects one has to shift from the lattice to a continuum 
picture. In a magnetic model with the spin configuration varying slowly over the lattice spacing, one can 
approximate using a long wavelength description, the d dimensional lattice to a continuum space R'^. At 
the same time the associated spin Sj would go to a vector field n"(x), a = 1, 2, 3 of unit length |np = 1, at 
the classical limit, after a proper renormalisatiou'^'' , which might induce topological invariants of different 
nature, at different dimensions d, under suitable boundary conditions 20-23 _ 

We intend to use this construction on a two-dimensional xy-plane perpendicular to the applied mag- 
netic field, to simulate the result on MnSi found through the topological quantum Hall experiment. The 
physically motivated boundary condition demands that at large distances: |x| = p ^ 00, the spin field 
n"(x) should go to its vacuum solution, orienting itself to a fixed vector = 6a3, along the applied 
field. This condition in turn introduces a nontrivial topology by identifying the infinities of the coordi- 
nate space to a single point, which compactifies the vector space to a sphere and defines thus a 
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mapping: —>■ S'^, linked to the nontrivial homotopy group 7r2(5^) = Z (see Fig. 1). Each topological 
sector is labeled by an integer A^-valued topological charge Q = J (i^x<I>(x), which defined as the degree 
of this mapping, shows how many times such a skyrmionic field configuration n(x) winds up or covers 
the target sphere, when the coordinate space is swept once. The corresponding charge density can be 
expressed through the spin vector field as 

cl>(x) = ^(n[5,nxa,n]). (1) 




Figure 1 | The induction of nontrivial topology. The two-dimensional coordinate space x G i?'^, 
due to the boundary condition limp^oo n = = (0, 0, 1), compactifies to a 2-sphere S'^, with its north- 
pole P identified to all points at the space- infinities. The unit vector field n(x) G S^, therefore describes 
a sphere to sphere mapping with the integer-valued topological charge Q = N, defined as the degree or 
the winding number of this mapping. 

Another important idea, that we use for constructing our magnetic model is to maintain a lower 
bound H > const.\Q\ for energy H of the system through its topological charge Q. Note that, such 
an energy bound would guarantee the stability of the finite energy solution in each topological sector, 
since the topological charge is conserved independent of the dynamics of the system and can not be 
changed by unwinding the spin configuration by smooth transformations. Moreover, when the energy 
reaches its lowest value saturating this bound: H = Hmin = const.\Q\, called Bogomolny limit'^'^, some 
extraordinary thing happens in almost all known models 17,18,21,23-25. j|. yjg^^ig g^j-^ exact solution to the 
field model, in some cases for all topological sectors. 

The unusual properties of the helical magnets like MnSi and especially their partial helical order 
created by the skyrmion spin states are believed to be due to the competing forces between the ferro- 
magnetic exchange and an effective Dzyaloshinskii-Moriya (DM) type interaction. The ferromagnetic 
interaction tries to align the spins parallel to each other, while the DM-term with the broken inversion 
symmetry tends to orient them perpendicular to each other, settling finally to a partial helical order 
of topological origin. Moreover, the skyrmionic spin texture in MnSi, as observed in the recent topo- 
logical Hall experiments, appears in the plane perpendicular to the applied magnetic field. Therefore 
to model such a system we construct our Hamiltonian H from a standard Heisenberg ferromagnet in 
two-dimensions: Hq = — '/X]<j,j'> ^jSj/, where the spin Sj , located at site j = (i, j) interact with its 
nearest-neighbour at j' = (i, j it l).(i ib 1, j), in the 2-dimensional lattice. We include next the DM term: 
Hdm = Ji X]<j j'>(-^j,j' ■ [^j ^ ^j'])' which in spite of its traditional form, bears some crucial differences in 
the structure of its Dzyaloshinskii (D)-vector Djj/, which would result to our field model allowing exact 
skyrmion solution. Firstly, we take the D-vector aligned along the applied field, which leaves nontrivial 
only its third-component Dfy, breaking the rotational symmetry in the internal spin-space. Our next 
deciding assumption is the nonlinear nature of this vector, given through a spin-dependent structure: 
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D^y = D^{S?,Sy) = f{S?){S? — Sy). Notc that the second factor obeys the usual antisymmetric ex- 
change of the D-vector, while the first factor, which we take in the explicit form f{S?) = (1 — {Sf)'^)~^ 
manifestly breaks the space-inversion symmetry of the vector: Dy^ ^ ^fi'' inducing the same property 
to the Hamiltonian of the system as usual through the DM interaction. 

Our total Hamiltonian, which has the required broken inversion-symmetry and the competing spin 
interactions of opposite trends, should now be tested for its topological properties, for which we have to 
go for the continuum limit S(j.j) n{x,y). For constructing the corresponding field model from the 2- 
dimensional lattice model, we have to take carefully the vanishing limit of the lattice constants in both the 
lattice directions which would result, for example, S(j-i-ij) — S(jj) — > ±dxn(^x,y)j ^{i,j±i) ~ ~^ 

±dyn(^x,y) 6tc. This procedure reduces our lattice spin-model H = Hq + Hdm to a field Hamiltonian 
with 

Ho = J J d'^x{Vnf and Hdm = Ji j d^x f{n^)Vn^ • (n^Vn^ - n^Vn^). (2) 

It is not difficult to see that, at this continuum limit the chirality of three non-coplanar spins x = 
Sj • [Sj' X Sj"] gets linked to the topological charge density ^{x), expressed through the n-field as ([T]). 

It is important to note, that by tuning the coupling constants J ^ Ji, our field Hamiltonian H, 
as can be shown by using a simple school-level geometric inequality^^ , would acquire the crucial lower 
bound : H > ATr^/3J\Q\, through the topological charge Q, the significance of which for the topological 
stability of the solutions is emphasized above. Focusing now on the conjecture, that the Bogomolny 
limit, at which the energy bound is saturated, should yield the exact finite-energy soliton solution, we 
look for the situation when this limit could be reached for nontrivial contributions from both parts of 
the Hamiltonian. We find that this is possible for a field configuration n*, where the space-inversion as 
well as the circular symmetry is lost. More precisely, the spin texture n* = (sin/3cos7,sin/3sin7,cos/3), 
should be such that, the directions V/? and V7, defining the helical order, would cross each other at an 
angle 60", with three field vectors: a = V/? , b = sin/?V7 and c = a — b, forming an isosceles triangle. 
Under this special geometric condition, the energy bound is exactly saturated attaining its lowest value 
-ffmm[n*] = 4Jtt\/^N, in each sector with Q = ±N. Analysing this intriguing geometry^^ we can extract 
the exact magnetic field solution as 

sin/? = 2(ppo)^et^/((p)^^ + (/>o)'^^e"^) and , 7 = ±Na, (3) 

in polar coordinates (/), a) with a constant scaling parameter pQ. This exact solution, which can be checked 
by direct insertion, obeys the required boundary condition n(p 00) = n^o = (0,0, 1), and yields the 
topological charge Q = zizN. The distribution of the charge density a), which coincides also with the 
energy density at the solution manifold n* is shown in Fig. 2, together with the magnetization sin /?(/?, a) 
and the spin component = cos I3{p, a) , calculated for the exact result specific to the first excited 
state with Q = 1. The figure shows that at the skyrmion axis p = 0, the charge density becomes maximum, 
vanishing gradually with the distance, while the magnetisation and the spin component exhibit a toroidal 
structure without circular and space- inversion symmetry. Their maximum/minimum is attained at a 

a 

distance pm = Po , varying with a, which describes an intriguing form of the torus, with its radius 
spreading out with increasing polar angle. As revealed in quantum Hall experiments^'^, the peculiarly 
large Hall conductance is created by an effective magnetic field proportional to the topological charge 
density. Therefore the detailed and exact form of the charge density together with the magnetization 
and the spin field component, found here based on the exact skyrmion solution, should serve as a guiding 
resource for the future precision experiment. 
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b c 
Figure 2 | Exact skyrmion state with Q = 1. a, The charge density ^{p,a) = s\v? [3{p,a) / . 
which coincides also with the energy density for the exact skyrmion. b, The magnetization |m| = 
((n^)^ + (n^)^)2 = sin /?(/), a), and c, Spin component = cos/3(p, a), showing toroidal structure with 
explicitly broken inversion and circular symmetry. 

The graphical result based on the numerical calculation on a model for MnSi, proposed in Ref 13, 
shows that the skyrmion lattice derived as a spontaneous ground state solution, preserves the circular 
as well as the space-inversion symmetry, when projected on two-dimensions. However, our exact result 
([3]) based on our magnetic model with the broken space-inversion and the spin rotational symmetry, 
induced by the DM interaction, gives a different picture. It shows, as depicted in Fig. 3, that the partial 
helical order of the spin texture on the xy-plane manifestly breaks both the circular and the inversion 
symmetry. The future experiments capable of identifying the finer details of the topological spin texture 
should be able to settle this issue. At the theoretical level, it is significant to note also that, the exact 
topological soliton without circular symmetry, presented here, is a nontrivial generalisation over the well 
known solution of the nonlinear cr-model^^'^^'^^, exhibiting this symmetry. In fact all other known exact 
solutions of the field models'^^"^^, obtained at the Bogomolny limit, are basically spherical symmetric, 
except only the present one. 

Another valuable information, e.g. the detail spin structure in higher topological sectors with arbitrary 
integer charge Q = N > 1, can be found easily from our exact solution Although the solution with 
higher topological charge was investigated by using a combination of analytic and numerical methods in 
connection with the Hall skyrmions ^ and recently for the bilayer graphine^^, no such study, as far as we 



know, is undertaken for the magnetic models. Therefore the present result, designed specifically for the 
magnets like MnSi. and more significantly having purely analytic result with manifestly broken inversion 
symmetry, is novel in many respect. Such an exact skyrmion with higher topological charge Q = 3, is 
shown in Fig. 3b. These stable states in higher sectors with more intricate texture, as we see in the 
figure, would be challenging to detect in topological Hall experiments. 

The topological soliton with non-circular symmetry, presented here, is exceptional as an exact field 
theoretic solution and would be interesting for its detection and applications in other systems, lacking 
such a symmetry. 




a b 
Figure 3 | Skyrmion spin texture with partial helical order, described by the exact solution 
of the spin vector field n(x). a , with the topological charge Q = 1 and b, Q = 3. Both the 
figures demonstrate manifestly broken space-inversion and the circular symmetry on the xy-plane, with 
a distinct difference in their helical pattern, describing different degrees of mapping in these two cases. 

References — 

[1] Bouick, M. J., Chandar, L., Schiff, E. A. & Srivastava, A. M. Cosmological Kibble mechanism in 
the laboratory: string formation in liquid crystal. Science 263 943-945 (1994) 

[2] Travers, A. DNA topology: dynamic DNA looping. Cwrreni Biology 16 R838-R840 (2006) 

[3] Bauerie, C. et al. Laboratory simulation of cosmic string formation in the early universe using 
superfluid ^ He. Nature 382 332-334 (1996) 

[4] Sondhi, S. L., Karlhede, A. &; Kivelson, S. A. Skyrmion and the crossover from integer to fractional 
quantum Hall effect at small Zeeman energies. Phys. Rev. B 47, 16419-19426 (1993). 

[5] Lillihook, D., Lejnell, K., Karlhede, A. & Sondhi, S. L. Quantum Hall skyrmions with higher topo- 
logical charge, cond-mat.9704.121 

[6] Al'Khawaja, U. k, Stoof H. T. C. Skyrmions in a ferromagnetic Bose-Einstein condensate. Nature 
411, 918-920 (2001) 

[7] Lee, M., Kang, W., Onose, Y., Tokura, Y.&; Ong, N. P. Unusual Hall effect anomaly in MnSi under 
pressure. Phys. Rev. Lett. 102, 186601 1-4 (2009). 



5 



[8] Neubauer, A. et. al. Topological Hall effect in the A phase of MnSi. Phys. Rev. Lett. 102, 186602 
1-4 (2009). 

[9] Pfleiderer, C, Julian, S. R. & Lonzarich, G. G. Non- Fermi-liquid nature of the normal state of 
itinerant-electron ferromagnets. Nature 414, 427 - 430 (2001). 

[10] Pfleiderer, C. Non-Fermi-liquid puzzle of MnSi. Physica B 328, 100 - 104 (2003). 

[11] Dolron-Leyraud, N. et. al. Fermi-liquid breakdown in the parametric phase of a pure metal. Nature 
425, 595 - 599 (2003). 

[12] Pfleiderer, C. et. al. Partial order in the non-Fermi-liquid phase of MnSi. Nature 427, 227 - 231 
(2004). 

[13] Ro/31er, U. K., Bogdanov, A. N. & Pfleiderer, C. Spontaneous skyrmion ground states in magnetic 
metals. Nature 442, 797 - 801 (2006). 

[14] Binz, B. & Vishwanath, A. Chirality induced anomalous - Hall effect in helical spin crystals. Physica 
5 403, 1336 - 1340 (2008). 

[15] Miihlbauer, S. et. al. Skyrmion lattice in a chiral magnet. Science 323, 915 - 919 (2009). 

[16] Skyrme, T. H. A nonlinear field theory, Proc. R. Soc. Lond. Ser. A 260, 127 - 138 (1961). 

[17] Belavin, A. A. & Polyakov, A. M. Metastable states of two-dimensional Heisenberg ferromagnet. 
JETP Lett. 22, 245 - 246 (1975). 

[18] Honerkamp, J., Patani, A., Schlindwein, M. & Shafi, Q. On kink solution in two space dimensions. 
Nuovo Cimento Lett. 15, 97-100 (1976). 

[19] Kundu, A. Instanton solutions in the anisotropic a - model. Phys. Lett. B 110, 61-63 (1982). 

[20] Fogedby, H.C. Classical and quantum aspects of the continuous Heisenberg chain at T=0. Z. Phys. 
5 41, 115-127 (1981) 

[21] Duff, M. J. & Isham, C. J. Soliton and vortex type solutions in nonlinear chiral theories. Nucl. Phys. 
B 108, 130 - 140 (1976). 

[22] Faddeev, L. & Niemi, A. J. Stable knot-like structures in classical field theory Nature 387, 58-61 
(1997) 

[23] Felsager, B & Leinaas, J. M. A conformally invariant cj-model in 4-dimensions. Phys. Lett. B 94, 
192-194 (1980) 

[24] Bogomolny, E. B. The stability of classical solutions. Sou. J. Nucl. Phys. 24, 449-454 (1976) 

[25] Actor, A. Rev. Mod. Phys. 51, 451-525 (1979) Classical solutions of SU{2) Yang-Mills theories. 

[26] Abanin, D. A., Parameswaran, S. A. & Sondhi, S. L. Charge 2e skyrmions in Bilayer graphine. 
.arXiv:0904.004 [cond-mat.mes-hall] (1 Apr. 2009) 



6 



